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I. INTRODUCTION 

In recent years modified dispersion relations which 
break Lorentz invariance have appeared in different con- 
texts of gravitational physics. 

On the one hand, they appeared in several works which 
address the so-called trans-Planckian problem [1-3]. In 
black hole physics the modes responsible for Hawking 
radiation reach arbitrarily high energies near the black 
hole horizon when measured by a free- falling observer. 
This observation has lead some authors to study the ro- 
bustness of Hawking radiation when modifying (from the 
outset) the dispersion relation beyond the Planck scale 
[4-7]. A similar problem appears in inflationary cosmol- 
ogy: the modes responsible for the large scale structure 
had length scales much smaller than the Planck length 
in the early stages of inflation. Therefore it is of in- 
terest to determine to what extent the properties of the 
fluctuation spectrum are sensitive to modiflcations of the 
dispersion relation at the Planck scale [8-11]. In these 
works, the introduction of a nontrivial dispersion rela- 
tion has been originally suggested by formal analogies 
with condensed matter [4, 12-15]. However, the main 
physical motivation is that unknown effects of quantum 
gravity (such as radiative corrections) might introduce 
nontrivial dispersion relations in these contexts. Follow- 
ing 't Hooft's [16] observation that strong gravitational 
interactions in the near horizon region might alter the 
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semiclassical description of black hole evaporation, a flrst 
dynamical realisation of this line of sought has been pur- 
sued in Refs. [17, 18]. 

On the other hand, there have been hints from high en- 
ergy astroparticle physics that we could already see the 
effects of a violation of the Lorentz symmetry. In particu- 
lar modified dispersion relations have been used in order 
to provide an explanation for observations of ultrahigh 
energetic cosmic rays (with energies higher than 10^" eV) 
beyond the Greisen-Zatsepin-Kuzmin (GZK) cutoff [19- 
22]. As in the trans-Planckian problem, several disper- 
sion relations have been exploited, inspired by various 
approaches to quantum gravity [23, 24], string theory 
[25], non-conmutative field theory [26], variation of cou- 
plings [27] and multiverses [28]. Modified dispersion re- 
lations often break the energy degeneracy for a given 
three-momentum: the particle energy can become helic- 
ity dependent [23, 24] and light can become birefringent 
[29, 30]. A popular form for these modified dispersion 
relations below the Planck scale is, in the cosmological 
rest frame, 

E'=m^ + \p\' + j:^r.^, (1) 

n>3 ^"Pl 

where Mpj is Planck's mass and rj^ are coefficients of or- 
der 1 if the Planck mass is the relevant scale. Jacobson et 
al. have studied the constraints on the possible values of 
the parameters 77„ based on current astrophysical obser- 
vational data [31-34]. It is important to notice that both 
the real and imaginary parts 773 are already observation- 
ally constrained to be much smaller than 1. This shows 
that the "natural" assumption that the rjn should be of 
order 1 is perhaps too naive. The possibility of quadratic 
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modifications of the dispersion relation is also strongly 
constrained [29, 30]. 

As we learned from particle physics, Lorentz invariance 
is a key element in renormalized quantum field theory [35, 
36]. Nevertheless, it might be broken by quantum gravity 
at a fundamental level, or, alternatively, it may be broken 
in an effective way in nontrivial backgrounds only. Indeed 
in backgrounds which possess a preferred reference frame, 
radiative corrections to sclf-cncrgies might contain terms 
which effectively break the (local) Lorentz invariance — 
i.e., the Lorentz invariance in the tangent plane. Let us 
emphasize that this second possibility does not imply any 
kind of fundamental breaking of the Lorentz symmetry 
or any new physics. Rather by "effective breaking of the 
Lorentz invariance" we mean that radiative corrections to 
the dispersion relation may contain terms which depend 
on vector or tensor fields characterizing the background. 

In this paper we explore a model in the framework of 
this second alternative. Wc consider the propagation of 
a scalar particle immersed in a thermal bath of gravi- 
tons. Our aim is to determine how the inertia of the 
thermal bath affects the propagation of the scalar parti- 
cle through radiative corrections. To this end, we shall 
compute the thermal corrections to its self-energy. The 
real part of the self-energy gives the thermal mass shift 
and might introduce as well a nontrivial dispersion rela- 
tion [37]. The imaginary part instead gives the rate of 
approach to thermal equilibrium [38]. 

In fact it is well known that quantum effects evalu- 
ated in nontrivial backgrounds induce non Lorentz in- 
variant terms in the dispersion relation. For instance, 
the self-energy of a charged fermion immersed in a ther- 
mal bath of photons has been computed in Rcfs. [37, 39- 
41], and the effects of the thermal bath on the speed of 
light have been studied in Refs. [42-44]. Let us mention 
also that electromagnetic interaction in curved space- 
times also leads to modifications of the dispersion re- 
lations. The effect of QED vacuum polarization on the 
speed of light in nontrivial spacetimes has been studied in 
Refs. [45 49]. However, most calculations done so far do 
not concern gravitational interactions. Gravity-induced 
corrections should be dominant at energies approaching 
the Planck scale, hence relevant for the situations con- 
cerned with the trans-Planckian problem. Additionally 
notice that gravitational interactions are universal and 
not limited to charged particles. 

Two important points concerning our approach should 
be mentioned. On the one hand, since gravitational inter- 
actions are non-renormalizable, the system of the scalar 
field coupled to gravity should be conceived as an ef- 
fective field theory [36, 50]. Although the full theory 
is non-renormalizable, low-energy predictions which do 
not depend on the Planck scale behavior of gravity can 
be extracted [51-53]. On the other hand, since thermal 
corrections do not affect the ultraviolet properties of the 
theory, thermal field theory fits well with the effective 
theory approach. At a technical level, we shall incorpo- 
rate the thermal effects through the real-time descrip- 



tion of thermal field theory [54-56] , which can be seen as 
a particular application of the closed-time path (OTP) 
method in field theory [57-59]. 

We should also emphasize that wc shall compute only 
the lowest order corrections. Hence, wc deal with free 
gravitons; i.e., we neglect the backreaction of the scalar 
field on the metric perturbations. The modifications 
to the metric propagator can be computed both in the 
framework of stochastic gravity [60 69] or in the large- 
N limit of quantum gravity [70] (both approaches can 
be seen to be equivalent [71, 72]). A computation along 
these lines (but in the semiclassical approximation) has 
been carried out by Borgman and Ford [73] . They study 
the fluctuations in the focusing of a bundle of geodesies, 
propagating in a spacetime with metric perturbations in- 
duced by a thermal scalar field. To our knowledge, the 
higher order corrections to the self-energy have not yet 
been computed in the quantum field theoretical frame- 
work. 

Gravity-mediated modifications of the dispersion re- 
lation have been also considered by Burgess et al. [74]. 
Working in the context of brane-world scenarios, these 
authors compute the change to the dispersion relation 
of photons and fcrmions generated by interaction with 
four-dimensional effective gravitons. In that case the 
primary source of violation of the Lorentz symmetry 
is the modification of the dispersion relation of effec- 
tive four-dimensional gravitons induced by certain extra- 
dimensional configurations. Similarly to our case, this 
modification is communicated to brane-bound particles 
through gravitational interactions. 

The plan of the paper is the following. In Sec. II we 
introduce the action for the system considered — namely, 
a scalar field coupled to gravity. In Sec. Ill we compute 
perturbatively the leading correction to the self-energy 
of the scalar field at zero temperature. In Sec. IV we as- 
sume that the scalar field and the gravitons are in ther- 
mal equilibrium at a given temperature, and compute the 
corrections to the retarded self-energy. Finally in Sec. V 
we analyze the results and discuss the physical signifi- 
cance of the terms which break Lorcnz invariance. The 
Appendixes contain reference material as well as some 
technical details of the calculation. In Appendix A we 
present the Feynman rules for the scalar field coupled to 
linearized gravity. In Appendix B we give the results of 
the calculation of some integrals in dimensional regular- 
ization. In Appendix C we give a brief account of the 
real-time approach to thermal field theory. Appendix 
D is mostly technical and is devoted to the calculation 
of certain integrals that appear in the calculation of the 
thermal contribution to the self-energy. 

Throughout this paper we shall use a system of units 
with /i = c = /cb = 1. The signature of the metric will 

be (-,+,+,+). 
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II. SYSTEM 

We consider a minimally coupled real scalar field cf) of 

mass m propagating in a spacctinic characterized by a 
metric g^^. The action for the field is 



1 

+ -m 



2 J,2 



and the action for the metric is 



(2a) 



(2b) 



where R is the Ricci scalar, g is the determinant of the 

and K = \^32ttG = \/327r Lp\ is the gravita- 



mctric, and k 
tional coupling constant, with G being Newton's con- 
stant and Lpi being Planck's length. Assuming that the 
metric is a small perturbation of Minkowski spacctimc, 
g^i, = rjfiiy + Khfj_,^, the complete action S = S'^^g + Sg can 
be decomposed into the free scalar field, graviton and in- 
teraction actions as S = Sc/) + Sh + Sint- These actions, 
expanded in powers of k, are 



2 j,2 



(3a) 



dS 



di.hd^h^'-' + -d^hdi^h] +0{k), 



(3b) 
(3c) 

where h = /i''^, T^uy is the stress tensor of the scalar field. 



and 



(5) 

Indices are raised and lowered with the background met- 
ric Tj^ij. We have kept only the free terms in the action 
for the gravitons Sh because we shall only compute the 
lowest order corrections of the self-energy of the (j) field. 

To compute these corrections we need to introduce 
counterterms in order to cancel divergences. Since our 
system is non-renormalizable, it has to be understood as 
an effective field theory, a low-energy approximation of a 
more fundamental theory at the Planck scale [36, 51-53]. 
In order to compute to a given precision where E 

is the energy of the process, one has to introduce all pos- 
sible counterterms compatible with the symmetry whose 



coefficients are of order k" at most. In our case the most 
general action for the counterterms for the scalar field ac- 
tion up to order which is compatible with the Poincare 
symmetry is 



-'count — 



d X 



1 



ml - m^)(j)'^ 



+ -(Z-l)(a^<^a''</) + mV) (6) 



where mo = m-l-O(K^) is the bare mass, Z = 1 + 0{k'^) is 
the field renormalization parameter and Co = C + 0{1) 
is a bare four-derivative; coefficient. The finite coefficient 
C is a priori unknown and constitutes an external in- 
put of the theory. The value of C should be determined 
by experiments or by knowledge of the underlying more 
fundamental theory. On the other hand, we do not intro- 
duce counterterms to the interaction action because we 
shall not compute vertex corrections. 

Additionally the graviton action (3b) must be supple- 
mented with a gauge-fixing term. We will work with the 
harmonic gauge d„h^'^ — {l/2)d^h = 0, whose appropri- 
ate gauge-fixing action is [75, 76] 

5gf = - 1 d^x (d.h'^-' - ^d^hj (^d^h^x - Id^hj . 

(7) 

No Faddeev-Popov ghost fields are needed since we will 
not consider graviton self-interactions. 



III. ZERO TEMPERATURE 

The aim in this section is to compute the leading con- 
tribution to the vacuum self-energy of the scalar particle. 
At zero temperature the self-energy S*^="'(p^) is related 
to the Feynman propagator G^""' (p) through 



(jTp 



332 -Hrn2 -FS(^=o)(p2)- 



(8) 



We recall that the self-energy can be computed as the 
sum of all one-particle irreducible diagrams with ampu- 
tated external legs [35, 36, 77]. In order to regulate di- 
vergences appearing in the calculation we will use dimen- 
sional regularization [78, 79]. 

The two diagrams which may contribute to order 
are shown in Fig. 1. One must also take into account the 
contribution of the counterterms: 



(9) 



We first concentrate on the first diagram S^^^ "'(p^)- Ap- 
plying the Feynman rules described in Appendix A we 
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p — k 





FIG. 1: The two Feynman diagrams needed for the calculation 
of the self-energy: respectively, E(i)(p) and S(2)(p)- 



find, in d spacetime dimensions 



(10) 



{p — kY — ie k^ + ' 



wiiere e = 4 — d, V^vais and T^v{p,k) are given in Ap- 
pendix A, and /x is an arbitrary mass scale. Developing 
the products in the numerator we find 



- K^m^(l + e/4)I{p) + 0{e), 



(11) 



where the momentum integrals I{p), I nip) and I^v{p) 
are defined and computed in Appendix B. The result for 

E(Tr°'(p2) is 



(47r)2 



'Ti „ A m p „ , 9 



+ 2m^p^ 



2p2 



4 , 2 2\ , IP^+m^ 
m + m p I m 



— ze 



(12) 



where 1/e = 2/e — 7-|-ln47r. The second diagram 
^(2f"'(P^) ^ massless tadpole, and these are identically 
zero in dimensional regularization [78]. 

Once we include the countertorms and take the limit 
£ ^ the renormalized self-energy is found to be 



In order to relate the bare and the renormalized param- 
eters we have imposed the on shell renormalization con- 
dition 



-i(T=0) 



(-m^) = 0, 



(15) 



which fixes the position of the pole to be m^. The com- 
plete specification of the on shell renormalization scheme 
also involves fixing the residue of the pole, which im- 
plies an additional condition on the derivative of the self- 
energy [35, 36]. However this second condition is found to 
be singular due to the divergence of the logarithms when 
evaluated on shell. This kind of singularity is generic 
for theories with interacting massless particles because 
of the absence of well defined asymptotic regions in the 
presence of long-range forces. For a discussion of this 
point in the case of QED see Ref. [35]. Anyway in the 
paper we are just interested position of the pole, which 
is always well defined. 

Note that the only divergence arises from the field 
renormalization. In particular, the four-derivative coeffi- 
cient C does not get renormalized: there are no logarith- 
mic corrections to the p* term in the self-energy (13), 
and C coincides with its classical value Cq. Therefore 
it would have been consistent not to include the four- 
derivative counterterm and simply take C = from the 
beginning. 

The self-energy develops a negative imaginary part for 

— p2 > TO^, 



ImEC-,(p2) ^ ^(_^2 _ ^2)_^ + ^ + -V 



(16) 



which, according to the optical theorem [77], accounts 
for the probability of a scalar particle with momentum p 
to decay into an on shell scalar and an on shell graviton. 
Notice that on shell (p^ = — m^) the imaginary part of 
the self-energy vanishes, because there is no phase space 
available for the spontaneous emission of a graviton. 



m 
T 



4 , 2 2\ ,„ / + 



(47r)2 V2p2 

- (i^K + ™V)ln 
+ CK\p'' + m^f + 0{K^) 



In ( 1 + — ie 
mr 



le 



(13) 



where we have absorbed divergences in the following way: 



Z = l + ICK^rn^ - 
Co = C + 0(k4). 



t^rr? ( 1 



(47r)2 



(14a) 
0{k% (14b) 
(14c) 



IV. FINITE TEMPERATURE 

In this section we compute the thermal contribution to 
the self-energy working within the real-time approach to 
thermal field theory. In this approach, which is briefiy 
explained in Appendix C, the number of degrees of free- 
dom is doubled, and one has to consider four propagators 
organized in a 2 x 2 matrix Gab{p)- The self-energy also 
becomes a matrix Ti°-^{p). We will concentrate on the 
retarded propagator Gr(p) = G\\{p) — Gi2{p) for the 
following reasons: First, it is the one that exhibits sim- 
ple analytical properties at finite temperature (analytic- 
ity in the upper p° plane); second, as shown in Appendix 
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C, it is directly connected with the retarded self-energy 
Sr = Eii(p) + through 



Gk{p) 



(17) 



and, finally, the position of its poles have well defined 
interpretations in terms of energies and thermalization 
rates. Furthermore, the retarded propagator is the one 
that is naturally obtained from analytic continuation 
from the Euclidean propagator in the imaginary-time for- 
malism. For more details we refer to Appendix C and to 
Refs. [54-56]. A quantum mechanical model where this 
point will be studied in detail will be presented elsewhere 
[80]. 



J 



A. Self-energy: Real part 

Let us proceed to the calculation of the real part of 
Sr(p). Since ReSR(p) = ReE"(p) [see Eq. (C19)], we 
shall compute the real part of instead. We have 

to consider both diagrams sketched in Fig. 1, where now 
internal propagators are taken to be thermal and of type 
11: 

Til 



E^^(p) = (m^ - m') + {Z- l){p^ + m^) + k'CqP^ 



S^i)b)+S[2)(f)+0(«^). 



The first diagram, Sj^^(y)), is given by 



(18) 



1 



{p — kY — 
1 

^2 _|_ _ 



+ 27rm(|p"-r|)^((p-fc)^) 



27rm(|/fc°|)(5(fc2-m2) 



9iip^,k'^,P- k) 



r 



(19) 



where 



9i{p^,k'^,P- k) = ^T^,/(p, k)Taf}{k,p)V 



= k^p^ - 2{k ■ p)m^ - 2(1 + s/4)m'^, 



(20) 



the function n{E) is the Bose-Einstein distribution, 

1 



n{E) 



1 - e^/T ' 



(21) 



and we recall that e = 4 — d and that r^y(p, k) and ■p''^"'' 
are given in Appendix A. On the other hand, at finite 
temperature the tadpole diagram S^2)(p) longer van- 
ishes since its temperature-dependent part gives a finite 
contribution, 

5^(2)(P) = - y53(p') / ^27rS{k')n{\k^\), (22) 



where 



2i 



93 



ip^) = --V^,„0{p,p)V^-'"^ = lOm^ + 4p2, (23) 



with V^vaisiPi k) given in Appendix A. 
To order one has 



ReS"(p) = ReI]'--°'(/) - ^ [^^^^ + + ^(^'')] 



+ 0{k% 



(24) 



where S<'^=°'(p^) is the T = self-energy [see Eqs. (9) 
and (13)] and where the integrals A{p), B{p) and C(p^) 
are defined through 



A{p) = j 



i{\p^-k°\)5{{p-kf) 



X gi{p',k^p-k)F- 



/ ^^ilk^rnk'+m') 
xgr{p^,k'^,p-k)V- 



ip-kr' 

C{p'')=93ip') I -^n{\k°\)S{k% 



(25a) 



(25b) 



(25c) 



with "P" meaning principal value. We have set d = 
4 because thermal contributions are ultraviolet finite. 
Note that the same renormalization process which makes 
the T = self-energy finite also renormalizes the 
T > self-energy; there is no need to introduce addi- 
tional temperature-dependent counterterms. Hereafter 
we will concentrate in the on shell results, = Ep = 
^Jw? + |p|2, since these are the ones we will need after- 
wards for the calculation of the thermal mass and the 
modified dispersion relation. 

The integrals A(j)) and Cij?) take into account the ef- 
fect of the thermal gravitons on the scalar particle. In the 
on shell case it is possible to give an explicit expression 
for these integrals valid at any temperature (the details 
of the calculation can be found in Appendix D): 

A(£;p,p) = -im2T2, (26) 
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a-m'^) = \m^T^. (27) 

Notice that A(i?p,p) does not depend on the three- 
momentum p. This is a surprise since Eq. (25a) iw not 
manifestly Lorentz invariant. This simplification, which 
is also found in electrodynamics [37] , will have important 
consequences when computing the modified dispersion 
relation. 

The integral B{p) takes into account the effect of the 

thermal scalars in the heat bath. For temperatures well 
below the mass m there are almost no scalar particles 
in the bath, and this indeed shows up in an exponential 



suppression of B{Ep, p) at low temperatures, T <^m: 

(28) 

For temperatures of the order of the mass m, there are 
thermal scalar particles in the bath and they give a sig- 
nificant contribution to the self-energy. At high temper- 
atures T ^ m, the leading contribution to the integral is 
given by 



S(-Bp,p) 



'W? + |pp 



24|p| 



In 



2Vm2 + |pP-|p| 
2Vm2 + |p|2 + |p| ; ■ 8 



+ -m^T^, T > m. 



(29) 



The details of the calculation for both the low- and high-temperature limits are given in Appendix D. 

Summarizing, according to Eqs. (24) and (26)-(29) the real part of the on shell self-energy in the low- and high- 
temperature regimes is given by: 



ReER(i;p,p)« { 



6 



5rp3 



87r3 



'2 .-m/T {rn^+^\pl_ 

3m2 + 4|p|2 
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^2 + |p|2 + |p| 



2v'm2 + |p|2 - IpI 



T<Cm, 
, T> m. 



(30) 



Notice that from the high-temperature result we can deduce that for massless particles the on shell self-energy is 
exactly zero. 



B. Self-energy: Imaginary part 



get 



Now we want to compute ImER (p). Similarly to the 
previous subsection we could compute T,^^{p) and then 
make use of Eq. (C19); however, it is somewhat easier 
for us to compute T,^'^{p) and T,'^^{p) and make use of 
relation (C18) instead. The self-energy T,^'^{p) is given 

by 

X n(p° - A:°) sign(p" - fcO) ^^^^ 
X SUp - kf)6{k^ + m2)5i(p2, k\p ■ k), 

where we used the property 0{—p^) + n(|p°|) = 
sign{p^)n{p^). In a similar way we find for E2^: 



■[l + n(A:°)]sign(fc°) 



2 J (27r)2' 
X [1 + n(p° - k°)] sign(p" - fc") 
X 5{{p - kf)5{k^ + m2)5i(p2, k\p- k), 



(32) 



where now wc used 0(p°) + JT-dp^l) = sign(p'')[l + n{p'^)]. 
Thus from Eq. (C18) and the two previous equations we 



ImER(p) = --£>(p), 

where D{p) is defined by 

D{p) = J -0^F{p^,k^)g,{p\k\p.k) 

X 6{{p-kf)6{k^ +m^) 

with 

F{p°, k°) = sign(fc°) sign(p" - fc°) 

X [l + n{k°)+n{p° -k°)] , 

which can be developed to give 



(33) 



(34) 



(35) 



Fip'', k°) = I sinh ( ^) csch ( csch ( 1^°' 



2T 



2T 



V 2r 



(36) 

After manipulating the integral D{p) (the details can be 

found in Appendix D) the imaginary part of the self- 
energy can be expressed as the following phase-space in- 
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tegral: 



lmT,B.{p) 



K' ni'{iir + 2p^) 
327r|p| 



rQ2 

/ dQF{p°,Q) 
jQi 



■ (37) 



bitrary temperature. For simplicity, we restrict ourselves 
to the casep'' > |p|, but including both (p°)^ > TO'^ + |p|^ 
and < + |pp. In this situation the integral can 

be performed analytically to give 



Let us now evaluate the integral in Eq. (37) for an ar- 



327r|p| 



■In 



^^^^ ' 4T(pit + |p|) 



sinh 



p")^-|p| 

4T(pii ~ 



SmJH^ 4T(pO_|p|) jSmn^^ 4T{pO + |p|) J_ 



(38) 



Making use of the property 

ln(sinha:;) w |x| - ln2 - Z7r6'(-a;) - e^^l^l, \x\ > 1, 
we easily obtain the low-temperature approximation 

,2 / ^6 gj^4 



K 



167r \2p 



2 2 



0{ 



,-rn/T 



), 



(39) 



which is in agreement with the zero-temperature result 
of Eq. (16). At zero temperature the imaginary part 
of the self-energy gives a net decay rate, but at finite 
temperature it is related to the rate F at which a particle 
approaches thermal equilibrium [38]: 



■ImEt 



(40) 



In this expression Fd and Fc are, respectively, the anni- 
hilation and creation rates of particles with energy 
which take into account stimulated absorbtion and emis- 
sion by the thermal bath. An arbitrary ensemble of par- 
ticles with distribution function f{E, t) approaches ther- 
mal equilibrium through 



f{E,t) = 



+ c{E) e 



-T(E)t 



(41) 



where c{E) depends on the initial conditions. Hence the 
imaginary part of the self-energy corresponds to an ob- 
servable quantity. 

To order a real particle can neither emit nor absorb 
a real graviton, because these processes are kinematically 
forbidden. Hence, at this order, the imaginary part of the 
self-energy must vanish on shell. However in the on-shell 
limit p^ ^ Ep = -^/ipp+rn?, we obtain, from Eq. (38), 



ImSR(p) 



■ In ' 



.o^Bp 327r|p| I ^m2 + |p|2- IpI 



(42) 



which is nonzero if T 7^ 0. 

This non zero result for the on shell self-energy is an 
artifact of not having introduced an infrared regulariza- 
tion; see Refs. [41, 81]. To illustrate this point one can 
take the on-shell limit directly in Eq. (37): in this limit 
Qii Q2 0, so that the phase space of the integral van- 
ishes while the thermal function F{p^ , 0) diverges. Hence 
the non vanishing result we obtained is a direct conse- 
quence of the infrared divergence of the Bose-Einstein 
distribution. When regulating the infrared behavior 
for instance by giving a tiny mass to the graviton — no 
imaginary part is found at p^ = —w?. Additionally, had 
we worked in an arbitrary gauge we could have verified 
that the result of Eq. (42) is not even gauge invariant 
[41, 81] [86]. 



V. DISCUSSION OF THE RESULTS 

At zero temperature, the position of the pole of the 
propagator gives the energy of the state and hence defines 
the dispersion relation, according to the Kallen-Lehmann 
specral representation [35, 36, 77]. The position of the 
pole is found to be 



(pO)2=m2 + |p|2 + E'"=°>(-^') 



+ |p|^ (43) 



where the second equality is a consequence of the renor- 

malization condition (15). The dispersion relation is 
clearly Lorentz invariant, as expected. 

Similarly, at finite temperature it can be shown [37] 
that the "effective" dispersion relation of the particle, 
given by the location of the poles of the retarded propa- 
gator, determines the inertial properties of the particle. 
The location of the poles is given by 



|p|2 = 



ReER,(p",p) 



(44) 



The thermal mass is obtained by setting p = 0: 

= -|-ReER(mT,0). (45) 

In a Lorentz invariant situation one would simply have 
(p0)2 = jji"^ + |pp, but in general there can be additional 
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dependence on the three-momentum p on the right-hand 
side, 



{py = + IpI^ + ^(«, T, Tot, p). 



(46) 



The Lorcntz-brcaking additional term in the dispersion 
relation leads to modifications of the group velocity of 
the particles v = dp°/dp: 



^ pO 2p0 dp ■ 



(47) 



Let us find the explicit form of the thermal mass and 

Lorentz-brcaking terms both in the low and high tem- 
perature regimes. Equation (44) can be solved perturba- 
tively: 



P\2 



RcEr(Sp,p) + 0(k4 



(48) 



where we recall that Ep = ■^/m'^ + |pp. At low tem- 
peratures the modified dispersion relation, according to 
Eqs. (30) and (48), is approximately given by 



0\2 



to| + |p|2 



-m/T 



3m2-h4|p|2 



(49) 



where the leading contribution to the thermal mass is 

n o 1 



TO^ 



TO 



(50) 



In the high three-momentum limit the term which mod- 
ifies the dispersion relation becomes a constant and can 
be reabsorbed in the thermal mass. The group velocity 
is modified according to Eq. (47): 



^2^9/22.3/2 g-m/T 

\/2^(3to2 + 4|p|2)2 



(51) 



Notice that at low temperatures the Lorentz breaking 
term carries a Boltzmann factor e^™/^. This is due to 
the mentioned fact that the nontrivial momentum de- 
pendence comes from the thermal scalar particles whose 
abundance is exponentially suppressed at low tempera- 
tures. Analogously to what happens in electrodynamics 
[37], the effect of the graviton bath only shows up in the 
thermal mass. 

At high temperature, T ^ to, the modified dispersion 
relation is found to be 



,0\2 



2 , I |2 
TOy + |p| 





^2Ep + 


IPh 


-1 


JpI ^ 


^2Sp- 


\P\) 





and the thermal mass is 
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^ 48 



(52) 



(53) 



At high three-momentum the modification of the disper- 
sion relation can be also reabsorbed in the thermal mass. 
The group velocity is given by 



1 + 



1 



K^T^TO^ 



96 (4m2-F3|p|2)|p|2 



1 ^27^2^4 



96 



■In 



2Er. 



2E„ 



(54) 



At high temperature the terms which break the Lorentz 
symmetry are no longer exponentially suppressed. Notice 
also that Eq. (52) shows that there is no modification to 
the dispersion relation for masslcss scalars. 

The Lorentz breaking term is negative in the low- 
temperature case, which implies that the speed of propa- 
gation is lowered with respect to the standard relativistic 
case. In contrast, the Lorentz breaking corrections have 
a positive sign at high temperatures, so that the speed 
of propagation is increased. However, it is always lower 
than the speed of light, as can be seen by expanding 
Eq. (54) in the ultrarelativistic limit: 



p / TOM48 + K2T2(ln3-l)l 
V « — 1 ^- K 



96|p| 



+ 



• (55) 



We should emphasize that the above breaking of the 
Lorentz symmetry is only effective: if one applies a 
Lorentz transformation simultaneously to the particle 
and the bath, all results are Lorentz invariant. This can 
be seen can by introducing the unit vector l^^ which gives 
the four-velocity of the thermal bath. Then the energy 
and the three-momentum of the particle with respect to 
the bath are 



P 



+ Pu.P'^ ■ 



(56) 



In this article for obvious reasons of simplicity we worked 
in the bath rest frame where = (1, 0, 0, 0). 

For low values of the momentum the effective disper- 
sion relation can be expanded in powers of p. For in- 
stance in the low-temperature case, for momenta satisfy- 
ing |p| < TO, the effective dispersion relation (49) can be 
expanded as 



(pY«mT + |pf 



+ 



^2^5/2^3/2 e-WT |p|2 

\/27r3 3to2 

^2^5/2^3/2 g-m/T 2|p|4 

TIP 27n? 



(57) 



It is worth mentioning the difference between this expan- 
sion and that of Eq. (1). In a fundamental approach to 
the Lorentz symmetry breaking one expects each addi- 
tional power of momentum to be suppressed by increas- 
ing powers of the Planck mass, whereas in the effective 
breaking approach we have pursued there are several en- 
ergy scales. Hence there is more freedom in the possible 
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FIG. 2: Diagrams such as the one depicted on the left would 
contribute to order to the thermalization rate at finite 
temperature, and hence to the imaginary part of the self- 
energy. However all those one- vertex diagrams are identically 
zero because of kinematical reasons. The first non-vanishing 
contribution would bo duo to Compton-Iiko diagrams such as 
the one depicted on the right, and hence it would be a 
contribution. 



values of the suppression factor. For instance, in the low- 
temperature regime — i.e., in Eq. (57) — the |p|'' term is 
suppressed by (T/m)^/^e~'"/'^ < 1 in Planck units. The 
fact that in the latter approach there are more energy 
scales also explains why modifications to the |pp term 
are present in Eq. (57) while they are not included in 
Eq. (1). Notice also that the corrections that we found 
only contain even powers of the momentum, and thus are 
in agreement with the results of Ref. [82] . 

In the present universe the effects we discussed are 
completely negligible when applied to electrons or pro- 
tons, both because they are proportional to the Planck 
length square and because they are exponentially sup- 
pressed. To obtain relevant effects, one should consider 
Planckian temperatures. However, in this case perturba- 
tion theory around a flat-spacetime approximation prob- 
ably fails. Therefore the results obtained here should 
be considered an indication that effective violations of 
Lorentz invariance indeed occur, and that gravitational 
interactions cannot be neglected when exploring highly 
energetic regions of dispersion relations. 

Whereas the real part of the self-energy gives the 
change in energy of the particle and hence the dispersion 
relation, the imaginary part accounts for the dissipative 
effects and gives the thermalization rate. As we have al- 
ready commented, at the order the imaginary part of 
the self-energy is zero. We do not see the thermalization 
effect since all processes which could contribute to it, such 
as the one shown in Fig. 2 (on the left), are kinematically 
forbidden. In order to account for thermalization effects 
we should compute the imaginary part of the self-energy 
to order k**. At this order, the particle can exchange mo- 
mentiim with the thermal bath through processes such as 
the one to the right of Fig. 2, which would correspond to 
the Compton scattering in electrodynamics. Notice here 
a significant difference with respect to the vacuum case: 
while at zero temperature the self-energy is purely real to 
all orders of perturbation theory (because the particle is 
stable), at finite temperature it will acquire an imaginary 
contribution to order k**. 

Let us end with a short summary of the main points 
of the paper. We have illustrated with a particular ex- 
ample how local Lorentz symmetry is effectively violated 
because of the interactions with a nontrivial ensemble of 



metric fluctuations, even if Lorentz symmetry holds at a 
fundamental level. We have also shown how this effective 
violation can be addressed at low energies within stan- 
dard physics. The main quantitative result of the paper 
are Eqs. (49) and (52), which explicitly show the modifi- 
cations of the dispersion relation. As in the electromag- 
netic case, this effect is exponentially suppressed at low 
temperatures. Moreover, the modifications of the disper- 
sion relation are suppressed when the three-momentum of 
the massive particle (defined in the heat bath rest frame) 
is much larger than the temperature and the mass. This 
last result is somewhat unexpected since the gravitational 
coupling grows with the energy, unlike the electromag- 
netic case. Therefore, in spite of the derivative coupling, 
no violation of Lorenz invariance is found in the high- 
momentum limit, at least at one loop. Finally we have 
also shown that gravitational interactions do not gener- 
ate a four-derivative term in the scalar field action. 
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APPENDIX A: FEYNMAN RULES 

At zero temperature the free propagator for the scalar 
field is 



(Al) 



and the free propagator for the gravitons in the harmonic 
gauge is [51, 52, 76] 



(A2) 



where 



'PlJ.uaff = ^ (vf^aVi^/S + rinpTlua - -^^^VnuVaP^ (A3) 

in d spacetime dimensions. In four dimensions, Vnvafi = 

{VnaVi'lS + VtilSrii^oc - V/J.i^'llocf}) /2. 

The two-scalar-one-graviton vertex, shown in Fig. 3, 
is given by [52] 



T^u{p, k) = —\Pfj.k„ + Puk^ - ip ■ k)rif^^ - rn'^Vfj.i^]i (A4) 
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FIG. 3: On the left, the two-scalar-one-graviton vertex, de- 
noted by Tf_iy(;p). On the right, the two-scalar-two-graviton 
vertex, denoted by Vfj,uap{p)- 



and the two-scalar-two-graviton vertex, also shown in 
Fig. 3, is given by [52] 



At finite temperature the free propagators for the 
scalar field are given by Eq. (Cll) and the free 
propagators for the gravitons are {^ab)ij.va0{p) = 
'Pf^i,af3Gfjj {p)\m=o- Vertices type 1 are described by r^^ 
and and vertices type 2 are described by — t^i/ and 

— Vfii^a^- 



APPENDIX B: INTEGRALS IN DIMENSIONAL 
REGULARIZATION 

In Eq. (11) we need the following integrals: 



where 



+ 2 iVlJ.i'Icxffpa- + Val3liJ.upa)P^k'^ 



(A5) 



crk 



1 



{2Try [k'^ + m2 - ie][{p - kf - ie] 
d'^k k„ 



r r ^ e f 'i'k 



, (Bl) 
f, (B2) 



e][{p-ky 



(B3) 



The result in arbitrary d dimensions after series expan- 
sion in £ = 4 — d is 



/ 

Ini^ip) 



(47r)2 

2(47r)2 



e 

1+2 

£ 



^lnfl + 4 



2(47r)2 



m 

K,2 



m 



In 1 



i€j—\n 

„2 



P 



»2 ^4 ^2 
1- 1 1 U _i_ 

2i 6i 6p'^ 6 9 



\6p^ 



m 
V 



In 1-H 



Pfj.Pf 
(47r)2 



1 

3i 



m 
3^ 



In 1 



P 



m 



p^ + 



-In 



+ 0{e), 



+ 0(£), 



6 



In 



-I- m 



. , 1 , / p^ 
— ^e — - In ' 



3p4 



m2 13 
6p2 18 



(B4) 
(B5) 

(B6) 



0(£), 



where 



1 2 

- = 7 + ln47r. 

£ £ 



r 



APPENDIX C: REAL-TIME THERMAL FIELD 
THEORY 



Correlation functions 



For the benefit of those readers who might be unfa- 
miliar with the subject we give in this appendix a brief 
account of those aspects of real-time thermal field theory 
which are needed in the paper. For a complete introduc- 
tion to the thermal field theory we address the reader to 
Refs. [54-56]; see also Ref. [83] for an approach similar 
to ours. 



In thermal field theory the initial state for the fields is 

assumed to be a thermal state at temperature T = 1//3, 
which is characterized by a density matrix 



-I3H 



Tr(e-M)' 



(CI) 
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*f Ret 
C3 



Y C4 



FIG. 4: Integration contour in the complex-time plane used 
in the real-time approach to thermal field theory. The choice 
(7 = 0"'" makes the formalism analogous to the OTP approach 
to field theory. 



where H is the Hamiltonian operator of the system (for 
the purposes of this appendix we shall consider a scalar 
field (j)). The generating functional Zc\j] is defined as 



Zc\j]=Tr(^pTce 



) , (C2) 



where (l>{x) is the field operator in the Hciscnbcrg picture, 
C is a certain path in the complex t plane, Tc means or- 
dering along this path and j{x) is a classical external 
source. By functional differentiation of the generating 
functional with respect to (/), path-ordered correlation 
functions can be obtained. The generating functional 
may be computed in a path-integral representation as 



Zc\j] =y"^0 giJcdtJ d''yi{C[4>{x)]+j{x)4>{x)}^ 



(C3) 



where C is the Lagrangian density and the contour C is 
now restricted to begin at some initial real time U and 
to end at ti — i^. The boundary conditions for the path 
integral are x) = — x). Different elections for 
C lead to different approaches to thermal field theory: a 
straight line from ti to ti — i0 leads to the imaginary-time 
formalism, and the contour shown in Fig. 4 leads to the 
real-time formalism. By choosing a = our formalism 
will coincide with the closed time path (CTP) approach 
to non equilibrium field theory [57-59]. In fact many of 
the properties listed below are not limited to the ther- 
mal case but are valid in a more general non equilibrium 
situation. We shall try to specify in each case which 
properties are general and which ones are particular to 
the thermal case. 

If we are interested in real-time correlation functions, 
the path along C3 and C4 can be usually neglected once 
we properly take into account the boundary conditions 
of the path integral. If we define ^1,2 (i,x) = 0(i,x) and 
ji^2(i;x) = j{t,x) for t e Ci^2) and take ti — > —00 and 
tf — !■ CX3 the generating functional can be reexpresscd as 
Zc[j] = Z^iZ[ji, j2]i where Z34 represents the constant 



contribution from the segments C3 and C4 and Z[ji,j2] is 



X e 



-2 J d'^X {C[<l>2{'<')]+Mx)Mx)} 



(C4) 



Correlation functions are defined by second functional 
differentiation with respect to the external source j{x): 



Gab{x,x') 



1 



Z[0,0] 6j''{x)6f{x) 



(C5) 



j 1=32=0 



where a, 6 = 1, 2. In operator language correlation func- 
tions can be written as 



Gl2{x,x') = {(t){x)(j){x')), 

G2i{x,x') = {${x')4>{x)), 
G22ix,x') = {Tcf){x)(}){x')), 



(C6a) 
(C6b) 
(C6c) 
(C6d) 



where (T) T is the (anti-) time-ordering operator, and 
where the average means (• • •) = Tr (p - • • ). Lowercase 
roman indices are raised and lowered with the "metric" 
Cab = diag(l, — 1). From the above expressions it can be 
readily seen that not all correlation functions are inde- 
pendent. The following relations in Fourier space are a 
consequence of Eqs. (C6): 

Gii(p) - GUp). GMp) = G2i{-p), (C7a) 
Gii(p) + G22(P) = Gi2(p) + G2i(p), (C7b) 

where correlation functions in Fourier space are 

Gabip) = j d^x e-'f^ GabiX + A/2,X- A/2), 

where we have introduced the new variables A = x — x' 
and X = {x + x')/2 [87]. Additionally, in the case of an 
initial thermal state the following relation, a consequence 
of the Kubo-Martin-Schwinger (KMS) formula (equiva- 
lent to the fluctuation-dissipation theorem, in another 
context), is also verified: 

Gn{p) + G22{p) = e-f^p" GM + e^^" G2i(p). (C8) 

Thus, in thermal field theory the knowledge of just one 
correlation function determines all of them. 

An important correlation function is the retarded prop- 
agator, which is defined as 

GR(a;, x') = eix'' - a;'O)([0(ar), ^{x')]), (C9) 

and is related to the other correlation functions through 

GR(p) = Gn(p)-Gi2(p). (CIO) 

The retarded propagator has the remarkable property 

that has well defined analyticity properties at finite tem- 
perature, as opposite to most other propagators: it is 
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analytic in the upper half of the complex p° plane. Fur- 
thermore, the retarded propagator is the one that one 
naturally obtains from an analytic continuation of the 
Euclidian propagator in the imaginary-time formalism 
[84]. 



2. Perturbation theory 

Perturbation theory can be organized in a similar way 
as in the zero-temperature case, but taking into account 
that at finite temperature there are two kind of vertices 
(1 and 2) and four kind of propagators (11, 12, 21 and 
22) which link the two vertices. Vertices type 2 carry 
an additional minus sign with respect to vertices type 
1. When computing Feynman diagrams one has to sum 
over all possible internal vertices. For a real scalar field 
of mass m, free propagators are given by 



/ 



p2 _|_ ^2 _ jg 



+ 2Tr5{p'^ +m'^)n{\p°\ 



p"^ + m? + ie 
1 1^ 
1 1, 



(Cll) 



where n{E) is the Bose-Einstein distribution function: 

1 



i{E) = — 



(C12) 



Thermal contributions to Feynman diagrams are always 
finite in the ultraviolet regime because the Bose-Einstein 
function acts as a soft cutoff for momenta larger than the 
temperature T. The counterterms which renormalize the 
theory at zero temperature also renormalize the theory 
at finite temperature. Note also that the thermal part 
of the propagator, which breaks the Lorentz symmetry 
through an explicit dependence on p°, is always on shell. 

Instead of working with the four propagators Gab {p) in 
the thermal case one can reorganize perturbation theory 
in a way such that just retarded and advanced propaga- 
tors arc involved [85]. In this case one has to consider 
how Feynman rules are transformed when working with 
the retarded/advanced basis. 



irreducible diagrams with amputated external legs that 
begin and end with type a and type b vertices, respec- 
tively. 

In general the self-energy components verify the non- 
perturbative relations 

e"(p) = -(e22)*(p), j:'^{p) = j:^\-p) (ci4a) 



(C14b) 



which can be obtained from Eqs. (C7) and (C13). The 
following equation is just verified if the initial state is 
thermal: 

S"(p) + S22(p) = - e-f'p" Ei2(p) - e^P" i:^\p). (C15) 

Thus, all the components of the self-energy can be deter- 
mined from knowledge of just one of them. Combining 
relations (C14) we obtain 



ImE^^(rt = -[E^^(p) + E^^(p)]. 



(C16) 



This last equation can be directly obtained from the cut- 
ting rules at finite temperature. 

A particularly useful combination is the retarded self- 
energy, defined as Er(p) = E^-^(p) -|- E^^(p). It is related 
to the retarded propagator through 



Gr{p) = 



p^ + m? 



(C17) 



The above relation, which justifies the name of retarded 
self-energy for Er(p), can be demonstrated by expand- 
ing the matrix equation (C13) and using the relations 
(C7) and (C14). Similar relations hold for the advanced 
propagator Ga{p) = (p) and the advanced self-energy 
Ea(p) — ^^r(-P)- Notice that a diagonal relation such 
as Eq. (C17) can be found only for the retarded (or ad- 
vanced) propagator. 

According to Eq. (CI 6), the imaginary part of the re- 
tarded self-energy can be also expressed as 

ImER(p) = ^[E2i(p)-Ei2(p)]. (C18) 
For a thermal state Sr(p) is related to S^^(p) through: 
Ei^(p) = ReE"(p) +tanh ImE"(p). (C19) 



3. Self-energy 

In a thermal or, more generally, in a non equilibrium 
situation, the self-energy has a matricial structure and is 
implicitly defined through the equation 

Gabip) = G<f^{p) + G^^]{p)[-ii:''\p)]GM, (C13) 

where G'fjj (p) are the free propagators of the theory. The 

ab component of the self-energy can be computed, simi- 
larly to the vacuum case, as the sum of all one-particle 



APPENDIX D: COMPUTATION OF A{p), B{p), 
dp'') AND D{p) 

In this appendix we shall compute the integrals A{p), 
B{p), C{p^) and D{p) which appear in the calculation of 
the self-energy at finite temperature; see Eqs. (25a)-(25c) 
and (34). 

Let us start by computing the integral A{p), defined in 
Eq. (25a). The Dirac delta can be expanded as 

6{{p - kf) = 6{q') = ^ [5{-q° + |q|) + 5(gO + |q|)] , 
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where we have introduced the new variable q = p — k sions, with 9 being the angle between p and q, and in- 
and where p = (p°,p) and q = (g°,q). Introducing now tegrating with respect to q^ with the aid of the delta 
spherical coordinates (c^, 6) in the three spatial dimen- function we get 



where Q = |q|, -P = |p|, a; = cos^ 



g2{p'',0,-p''Q + PQx) 



+ P- 



52(p',0,pOg + PQar) 



-(pO)2 + p2 + 2pOQ-2PQa; + m2 -{p'^)'^ + P'^ - 2p^Q - 2PQx + m'^ 



32(p^,g^P • q) = (P - qf,P- (p-q)) = -2m^ - 2mV + / + 2m^(p • q) - 2p^{p ■ q) +p^q^, 

and we have performed the trivial angular integration over (j). We now integrate with respect to x to get 

n{Q)dQ 



A{p) 







167r2p 



%PQ[-m' - (p") + P ] + m'{m' + 2[{py - P']} 

[m^ - (p" - P + 2Q)(p° + P)] [m^ - (p" - P)(pO + P - 2Q)] 
[m2 -{po-p- 2Q){p° + P)] [m2 - {pO - p){pO + p + 2Q)] 

I 



X In 



The result of this integral cannot be given in closed 

analytic form, in general. However in this paper we 
are mainly interested in its on-shell value p° = Ep = 
V m? + P2 , and in this limit the integral can be com- 
puted exactly at any temperature since the logarithmic 
term in Eq. (D2) vanishes. In this case the value of the 
integral is given by 

A{Ep,^) = J dQn{Q)Q = --m^T^, (D3) 

Jo 



where we used that 



f 

Jo 



dQn{Q)Q = 



(Dl) 



(D2) 



We now proceed with the computation of P(p), defined 
in Eq. (25b). Repeating similar steps as in the previous 
integral we get 



B{P) = f 
Jo 



n{E^)K'^AK 
2{2ttYE^ 



gi{p^,i):-p°Ek + PKx) 
-ijfif + P2 + 2p°E]^ - 2PKx 

I 



gi{p^,Q,p^E]^ + PKx) 



-(p0)2 + P2 - 2pO£;k - 2PKx 



where we recall that i^k = ^/rn? + K"^. The integral with 
respect to x can be analytically performed; the result is 
a very large and cumbersome expression, which we shall 
not reproduce here. The resulting expression cannot be 
integrated again in a closed analytic form. However, we 
may find particular expressions valid at low and high 
temperatures. As in the case of previous integral, we will 
restrict to the on shell results. 

At low temperatures, only those momenta k whose cor- 
responding energies are at most of the order of the tem- 
perature, Pk ^ P, contribute significatively to the inte- 
gral because of the presence of the thermal factor n(Pk), 
which acts as a soft cutoff. Hence, low temperature also 
implies low-energy and low momentum. Therefore, in the 
low-temperature approximation we may retain only the 



leading term m a K expansion: 
m2 + 2P2 



B{Ep,^) 



7r2(3m2 + 4P2) 

pOO 

X / dKn{E]^) [K'^ + 0{K^) 
Jo 

Taking into account that for low temperatures 
n{E^) « e-™/^ e-KV(2mT) 

and that 

dK e-^V(2mT)^2 ^ y|(^r)3/2^ 
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we find the following expression for B{p) at low temper- 
ature: 



5y3 



2P2 



2n-'> \ 3m 2 + 4P2 



-m/T 



(D4) 



Since the leading contribution to the integral is given 
in the ultrarelativistic regime, we can approximate the 
energy by the momentum in the Bose-Einstein function, 
n(£'k) ~ n{K). With this approximation we find 



We have not made precise the exact meaning of the "low- 
temperature" approximation employed above. In prin- 
ciple this approximation would require the temperature 
T to be much smaller than any relevant quantity with 
dimensions of energy that could be formed by a combi- 
nation of m and P. However, a detailed analysis of the 
expressions shows that the condition m » T is sufficient 
to guarantee the validity of the result. 

Let us now proceed to the calculation of B[p) in the 
high temperature regime. Since B{p) would be divergent 
if no thermal cutoff were present, at high temperatures 
the leading contribution to the integral is given by those 
momenta close to the temperature T.Thus as a first ap- 
proximation we can retain only the leading term mal/K 
expansion: 



B{Ep,p)= „.„P ln' P 



24P 



2Ep + P 



+lm'T'. (D6) 



Analogously to the low-temperature case, the high tem- 
perature approximation would a priori require the tem- 
perature T to be much higher than m, P and any relevant 
energy scale formed by combination of these two. Again, 
one can show that the condition T » m is sufficient to 
guarantee the validity of Eq. (D6) . 

We now move to the integral C{p^), defined in 
Eq. (25c). Its evaluation is straightforward: 



C{p') = gsip') / 
Jo 



KdK T"^ ,,r. 2 . 2X 

-^n{K) = —{lQm^ + 4.p^). 



B{E^,p) 



47r2p 



In 



2\frr? 



2Vm2 + P^ + P 



POO 

/ dKn{Ek)[K + 0{l/K'') 
Jo 



(D7) 

We only need its on shell value C{—m?) ~ T'^rr? 1%. 

Let us now consider the integral D{p), defined in 
Eq. (34). We start by introducing the variable q = p— k: 



Dip) 



I 



dg° d^q 
(27r)2 



X 5{ - {qy + q2) 5{ - (p" - q^ + (p - q)^ + m% 
where p= p) and q= {q^,c[). Next we expand the ffist delta function and integrate over q^: 



+ 



f d^q 
J (27r)3 



27r 



(D8) 



2|q| 



-|q|) 52(|pP - (/)^ 0, p • q + /|q|) 6{ - (/ + |q|)2 + (p - q)2 + m^). 



We now introduce spherical coordinates over q and expand the second delta function: 

Dip) = —j^dxJ^ dQ P(p°, Q) g2(P^ - (/)^ 0, QPx - Qp'^) 6{x - i/ + 2pOQ + m')/i2PQ)) 

+ ^ y ^ da; dQ Fip°, -Q) g^iP^ - ip^, 0, QPx + Qp'') 5{x - ip^ - 2p°Q + m^)/i2PQ)). 



where Q = |q|, P = |p| and x = cos 6*. We have al- 
ready performed the trivial angular integration over (f). 
Integrating with respect to x with the aid of the delta 
function we get 



Dip) = 



fl2(p^0,(p^+m2)/2) 
SttP 



rQ2 

/ dQFip°,Q) 



(D9) 



with 



Q2 = 



,0^2 



P2 



2(pO - P) ' 
which can be finally arranged as 
-m2(m2 -I- 2p2 



Qi = 



P\2 



2(p0 + P) 



Dip) = 



SttP 



Qi 



dQFip^Q) 



(DIO) 
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Recall that = -{jpf + P^. 
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